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Stellar population gradients from cosmological simulations: 
dependence on mass and environment in local galaxies 
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ABSTRACT 

The age and metallicity gradients for a sample of group and cluster galaxies from N- 
body-t-hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf 
galaxies show null age gradient with a tail of high and positive values for systems in 
groups and cluster outskirts. Massive systems have generally zero age gradients which 
turn to positive for the most massive ones. Metallicity gradients are distributed around 
zero in dwarf galaxies and become more negative with mass; massive galaxies have 
steeper negative metallicity gradients, but the trend flatten with mass. In particular, 
fossil groups are characterized by a tighter distribution of both age and metallicity 
gradients. We find a good agreement with both local observations and independent 
simulations. The results are also discussed in terms of the central age and metallicity, 
as well as the total colour, specific star formation and velocity dispersion. 

Key words: galaxies : evolution - galaxies : galaxies : general - galaxies : elliptical 
and lenticular, cD. 



1 INTRODUCTION 

Radial profiles of colours, ages and metallicities of stel- 
lar populations have been shown to be efficient tools to 
discriminate among different galaxy formation scenarios 
across a wide range of masses and environments. Ear- 
lier observational works were unable to assess the depen- 
dence of gradients on mass, mainly becaus e of the lim- 
ited s amples studied dPeletier et alj 1990al. Davies et alj 



iTamura fc OhteJl200d. l2003t). although it might have a non 
marginal role (e.g. lLa Barbera et al. 120051 ). 



19931, iKobavashi fc Arimotol Il999t ITamura fc Qhtal 120031 
La Barbera et alj 120051 ). Only recently, trends of gradients 



with mass have been more firmly assessed (e.g. iForbes et al.l 
2005), pointing towards different physical mechanisms forg- 
ing the behaviour of more and less massive g alaxies 



|Spolaor et alj|2009l . iRawle et al.1l2010l . iTortora et afll20ld ) . 

However, the dependence of the stellar population gradients 
on environment is still controversial l|Tamura et al ] l2000l ; 
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Simulations of galaxy formation are the ideal tool to 
interpret observations in terms of the underlying physi- 
cal mechanisms driving the evolutive processes of galaxies 
in various environments and a large mass range. For in- 
stance, in ITortora et alj (|2010h (T+10 hereafter) we have 
compared the stellar population (age and metallicity) gra- 
dients of a wide sample of local galaxies extracted from 
the Sloan Digital Sky Survey, with a number of literature 
simulations (of both monolithically collapsing systems and 
merger remnants). We have confirmed the metallicity as 
the main factor in shaping the colour gradient as a func- 
tion of mass. However, age gradients are still important, 
particularly for younger galaxies. We have shown that less 
massive systems, which are unlikely to get assembled by 
merging, have experienced a simple monolithic collapse, re- 
sulting into a nega t ive and steep m etallicity gradient (e.g. 
iLarsoni Il974 1 19751 . ICarlberd 1 1984 ) which decreases with 
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Figure 1. Age (left panels) and metallicity (right panels) gradients as a function of stellar mass (assuming an AY IMF). Top panels. 
Black filled and black open circles are for cluster and group members. The lines connected by gray triangles and green boxes are the 
medians for groups and clusters, respectively. The error bars are the quartilcs of sample distribution in each mass bin. At the bottom 
of the panels we also show the uncertainties on galaxies at different stellar masses, shown as error bars and determined as described 
in the text. On the right panel we also show the results from local cluster (green thick lines) and groups (gray thick lines) galaxies in 
ISpolaor et al ] (|2009h . the polygonals bound the regions where these data are distributed. Bottom panels. Black open boxes, gray open 
stars, red and blue points are for galaxies in normal and fossil groups, inner and external regions of clusters, respectiv ely. Our results are 
compa red with findings from other simulations. Th e pred icted metallicity gradients from dissipative collapse model in iKawata fc Gibson! 
( r 2003") and those from mergings in iBekki fc Shioval (119991") are shown as lo ng- and short- dashed lines. The shaded region is for remnants 
of major-mcrgings between gas-rich disk galaxies in lHopkins et al" (2009a). Blue and pink lines are the results from the chemo-dynamical 
model in Kawata (200l|), respectively for strong and weak SN feedback. 



mass (jKawata fc Gibsonl lioO&l IPipino et alj|20ich . Super- 
novae (SN) feedback is strong in this mass regime, reaching 
its maximum efficiency in the lowest mass syst ems and this 
contr ibutes to produce such negative gradients |Pipino et al.l 
2008). On the higher mass side instead, the trend with mass 
is inverted, demanding the contribut ion of strong joint ef- 
fect of mergings and AGN feedback dBekki fc Shioya 1999, 
iKobavashil 12004 iHopkins et"aTll2009al. ISiiacki et all 1 20071) . 
In spite of these results, recently IPipino et al .1 (|2010l ~ have 
reproduced the flattening of metallicity gradients in massive 
systems within a monolithic scenario of galaxy formation, 
pointing to a fundamental role of star formation efficiency 
in shaping the trend with mass. 

In the present paper we will use the TreeSPH 
N-body+hyd ro dynam ical simulations describ ed in 
iRomeo et all lj200rj l and iRomeo et all (120081 ) (here- 
after R+06 and R+08) to derive the profiles of both age 
and metallicity inside the galaxies and discuss them as a 



function of stellar mass and environment. This smoothed 
particle hydrodynamics (SPH) simulation has some advan- 
tages when compared with either semi-analytical models 
(SAM) or pure hydrodynamical simulations of individual 
objects. In general SAMs make use of large scale N-body 
simulations and rely upon different sets of assumptions and 
approximations (e.g. dark matter haloes are spherically 
symmetric and infalling gas is shock-heated to the virial 
temperature of the halo) that depend on many intercon- 
nected and variable parameters (such as star formation and 
feedback efficiency or IMF). In this scheme, the evolution 
of galaxies is adjusted on top of the pure collisionless DM 
haloes, by means of merger trees and models of synthetic 
colours. On the other hand SPH simulations contain fewer 
simplifying assumptions (as, for example, no restrictions 
on halo geometry), but have to restrict themselves on 
smaller scales (namely, galactic) in order to maintain a 
high resolution, resulting in a more limited dynamical 
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range. However, in terms of gas cooli ng both approaches 
are confirmed to give similar results ll Benson et al.l l200ll; 
iHellv et alJl2002l ). Our approach roughly lies in between, 
since the formation of galaxies is followed ab initio within 
a large cosmological volume, thus accounting for their 
interaction with a large scale environment. At the same 
time, it allows for describing in a self-consistent way the 
mutual cycle between inter-galactic medium, star formation 
and stellar feedback at the cluster scale, by means of 
the baryonic physics implemented in the hydrodynamical 
code, in particular the chemical enrichment of the gas 
surrounding galaxies (for details, see R+06). This is though 
attained at expenses of the final stellar resolution, that 
cannot reach the level of either pure N-body simulations 
used in SAMs, nor of the single-object hydrodynamical 
simulations at smaller scale. In this paper then we will test 
our SPH simulations against both observations and models, 
with the caveat of latter's different resolution: e.g. those o f 
gas- rich mergers between disks by iBekki fc Shioval || 19991 ) 
and lHopkins et all (|2009aT l. or the cosmo logical simulations 
includ ing a chemo-dynamical model by iKawata fc Gibson! 
l|2003h . 

In ij2]we will present the simulation setup and the fitting 
procedure adopted to recover the gradients; systematics in 
the fitting procedure and simulation resolution are analyzed 
in JjA] In §3] and [4] the trends with mass and environment 
are analyzed and we give an interpretation of the physical 
processes, while the conclusions are drawn in fj5] 



2 SIMULATION SETUP 

We have extracted simulated galaxies from the SPH sim- 
ulations in R+06 of two clusters of temperatures T ~ 3 
(CI) and 6keV (C2) and 12 groups (T ~ 1.5 keV), four of 
which are fossil. They were drawn and resimulated from a 
DM-only cosmological simulati on run with the code FLY 
l|Antonuccio-Delogu et alJlioolj ). for a standard flat ACDM 
cosmological model (h — 0.7, Q m = 0.3, as = 0.9) with 
150/i _1 Mpc boxlength. When resimulating with the hy- 
drocode, baryonic particles were added to the original DM 
ones, which were split according to a chosen baryon fraction 
ft = 0.12. 

Galaxies are composed by A r par star particles bound to 
the DM halo; each star particle represents a Single Stel- 
lar Population (SSP) of total stellar mass corresponding 
to the stellar mass resolution (M+^ssp) of the simulation. 
This is M^^ssp = 3.1 x 10 7 /i _1 for groups and CI, while 
M+ iSS p = 25 x 10 7 /i _1 for C2. The total luminosity and 
mass for each galaxy is defined as the sum of the luminosities 
and masses of their N par star particles. The individual stel- 
lar ma sses are distributed according to an lArimoto fc Yoshiil 
l| 19871 . hereafter AY) IMF; each of these SSPs is charac- 
terized by its age and metallicity (Z), from which lumi- 
nosities are computed by mass-weighted integration of the 
Padova isochrones. The "standard" super-wind model for 
SN is adopted, i.e. a prescription for SN feedback in which 
70% of the energy feedback from SN type II goes into driving 
galactic super- winds (see R+06 for more details). Although 
AGN feedback is not taken into account in this simulation, 
R+08 has shown that the colour properties of galaxies are 



fairly good reproduced; however, AGNs would play a major 
role only in high mass systems. 

Because of the galactic winds expelling baryons out into 
the IGM, lower mass galaxies have a higher fraction of DM 
over stars (and gas), which also results in a low absolute 
number of star particles. Throughout this paper we will 
only deal with the stellar component, willing to remain as 
much as consistent with the approach followed in the obser- 
vations. For this reason it is important to determine a lower 
limit to the number of star particles, in order to avoid sys- 
tems whose stellar component is not sufficiently resolved. We 
have originally applied a completeness limit in stellar mass 
of logM* > 9.5 Mq, but to be more conservative we have 
only retained those systems having N par > 80. We are left 
with a sample including 32 galaxies from the two clusters 
and 97 from the 12 groups, at z~0. The systematics in the 
galaxies with a lower number of particles will be discussed 
in § IA2I where we also show how the results could change 
when systems with N par ^ 80 were included. 

Brightest central galaxies (BCG) are excluded from this 
analysis, since their diffuse envelope mixing with the intra- 
cluster light makes difficult to fit a reliable luminosity profile 
on the basis of the galaxy-bound star particles only. 

For a given stellar population parameter X (= age, Z), 
we will assume the radial variation as logA(_R) = ax + 
bx log R, and define the X gradient as the angular coefficient 
bx = Vx = ^i'logfl^ • F° r convenience, we adopt as central 
quantity the value of log X at R — 1 kpc, i.e. the intercept 
ax- By definition, if the gradient is positive, i.e. Vx > 0, 
then X is lower in the central regions, while X decreases as 
a function of R when the gradient is negative, i.e. Vx < 0. 

The centre of the galaxy is defined as the median of 
the positions of all the star particles. The fit is made after 
having discarded those particles with age and metallicities 
> 3cr away from the median (see § Ali tor further details). To 
take into account disturbed galaxy shapes, we performed the 
log-log fit on three projected orthogonal planes and took the 
median of the parameters to obtain the final gradient values. 



3 RESULTS 

We show the age and metallicity gradients for groups and 
clusters as a function of stellar mass in Fig. [1] As shown in 
the upper panels of this figure, on average, low mass systems 
have null age gradients and very shallow (either negative or 
positive) metallicity gradients. However, a cloud of high age 
gradients in galaxies with log M* < 10.5 Mq is found in the 
external regions of simulated clusters and in normal groups, 
while galaxies in cluster cores and fossil groups have a quite 
tight distribution around V agc < 0.05 (see bottom panels). 
This is an interesting hint of the effect of the environment 
on the low mass systems, acting on flattening age gradients 
in higher density regions. Moreover, galaxies in the groups 
have shallower metallicity gradients, when compared with 
clusters. In particular, FG galaxies are distributed along a 
tighter sequence in both the age and the metallicity gradi- 
ents. 

At higher masses we have found slightly positive age 
gradients and negative metallicity gradients. On the side 
of the most massive galaxies (logM* :>1O.7M0), the mean 
age gradients become increasingly positive with mass at the 
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Figure 2. Age (top panels) and metallicity (bottom panels) gradients as a function of stellar mass. Filled and open symbols are for 
cluster and group galaxies, respectively. Continue gray and green lines are the medians for groups and clusters. From the left to right we 
divide galaxies in bins of total colour B — V , specific star formation SSFR, and velocity dispersion. First panels: Red and blue symbols 
are for galaxies redder and bluer than B — V = 0.8. Second panels: Red and blue symbols are for galaxies with log SSFR lower and 
higher than —2. Third panels: Red and blue symbols are for a total velocity dispersion higher and lower than 150km/s, respectively. 



high-mass end; on the contrary, the trend of metallicity gra- 
dients at high masses gets flatter with mass. At these mass 
scales, (normal) group galaxies show higher positive age gra- 
dients and steeper negative metallicity gradients. In particu- 
lar, while groups present a continuous steepening of metallic- 
ity gradients with mass, the mass trend of cluster members 
turns from steep to fiat at around logM* ~ 10.3 — 10.4, to 
slightly increasing thence on. 

Such trend inversion in the slope of metallicity gradi- 
ents is corroborated by comparing with data from a sample 
of cluster and group galaxies collected in lSpolaor et al.ll2009l 
(see also T+09) and shown in the top right panel of Fig. 
[T] Although their low mass region is populated by cluster 
galaxies only, at the massive side cluster galaxies have shal- 
lower metallicity gradients than the ones in group galaxies, 
qualitatively confirming our trends. The slope of the mas- 
sive side is even steeper with mass than ours, probably and 
partially due t o a dearth of very m assive galaxies in our sim- 
ulated sample IjRomeo et al.ll2005l ). apart of the BCGs. Like- 
wise, our cluster data follow a less steep trend than theirs 
at the low mass side too. 

As shown from the error bars in the upper panels of 
Fig. [T] the estimated gradients in low mass systems are in- 
trinsically noisier than the ones for massive galaxies. The 
error bars give an estimate of the uncertainty in the gradi- 
ents and contribute to the scatter in the derived trends with 
stellar mass. See i]A2l for further details about the procedure 
performed to obtain such results. 

In the same Fig. [T] we add a comparison with other 



models as well, mostly from hig h-resolution hydrodyn ami- 
cal simulations of di sk mergers (|Bekki fc Shiovalll999l and 
iHopkins et al.ll2009al ). In particular, the Z gradients of mas- 
sive systems fa l l in the range covered by merger remnants in 
I Hopkins et alj |2009al ). Our trend for metallicity gradients 
is also qualitatively consiste nt with the predic t ion o f the 
chemo-dynamical model by iKawata &: Gibson] (|2003r i. In- 
stead, our results are quite inconsistent wit h the almost flat 
gradie nts derived from merging model in iBekki fc Shioval 
(1999). Finally, when comparing with simulations in lKawatal 
(2001), we note that our points at intermediate mass are 
placed along their curve derived from models with strong 
stellar feedback, at least till 10 M©. To this regard, in the 
discussion below we will come back upon the role of SN 
winds in shaping the trend of Z gradients with mass. 

To check the dependence on other parameters than stel- 
lar mass, in Fig.[2]we discuss the trends in Fig.[T]by splitting 
the sample in two classes of galaxy colours, specific SFR over 
the last Gyr (SSFR = SFR/iW*) and velocity dispersion. At 
log M* <; 10.5 the bluest galaxies have the highest (posi- 
tive) metallicity gradients, while this result is inverted for 
massive galaxies. The age gradients do not show particular 
trends. In terms of SSFR, only few galaxies are star forming 
at z=0 (see R+08), and mostly the more massive ones in 
groups; these have positive age gradients and the steepest 
negative metallicity gradients. Note that few of the most 
massive galaxies (with logM, Jjl0.7) in our sample, with 
a recent star formation history (log SSFR > —2) and bluer 
colours (B — V < 0.8), present steep age and metallicity gra- 
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Figure 3. Age (top panels) and metallicity (bottom panels) gradients as a function of stellar mass. Filled and open symbols are for 
cluster and group galaxies, respectively. The continue gray and green lines are the medians for groups and clusters. From the left to right 
we divide galaxies in bins of central metallicity and age. Left panels: Red and blue symbols are for galaxies more and less metal-rich than 
0.008 at their centre. Right panels: Red and blue symbols are for galaxies centrally older and younger than 6Gyr. Here, we compare with 
results from the analysis of the local sample of ETGs in T+10. The boxes with bars are the results for all the sample, while the squares 
and triangles are for galaxies with a central age < 6 Gyr and > 6 Gyr, respectively. 



dients and are likely candidate for late-type systems. Merg- 
ing and SN feedback alone do not allow these systems to 
move onto the red-sequence, requiring the inclusion of other 
sources of feedback (e.g., AGN feedback) to quench SF and 
flatten the gradients. Finally, we also analyze the effect of a 
cut in velocity dispersion: due to the tight correlation with 
stellar mass, the systems affected by this criterion are pre- 
dominantly low mass systems in groups. 

In Fig. [31 galaxies are classified by their values of 
central metallicity and age. We confirm that both metal- 
licity and age gr adients depend on central quantities, re- 
spectively (e.g., iHopkins et"afl l2009al . IRawle et all l20ld 
T+10). In particular, metallicity gradients strongly depend 
on the central metallicity, with metal-poor systems hav- 
ing higher values (null or positive). In the last panels, we 
can also report a substantial agreement between the sim- 
ulations and the median gradients recovered from local 
ETGs in T+10 (black symbols in the right panels of Fig. 
[3]). As found in T+10, here we can see that metallicity 
and age gradients are strong functions of central age: mas- 
sive and older systems have null age gradients and shal- 
lower metallicity gradients (~ —0.2, —0.3). This is quite 
in agreement with other observations (e.g. ISpolaor et al.l 



120091 : IRawle et alJl201Ch a nd simulations of massive merger 
remnants discussed above ( Bckki fc Shiova|[l999l : lKobavashil 
120041 : IHopkins et alj|2009af ). On the other hand, the agree- 
ment between observed and simulated younger galaxies is 
looser, what might be tracked to the absence of those field 
galaxies in the sim ulated sample that are expected to have 
steepe r gradients (|La Barbera et all 120051 . iTortora et al.l 
l2010bl ). 



4 DISCUSSION 

The net result of the observed quantities at z~0 is the re- 
flection of all the physical processes modelled in the simula- 
tion code (namely radiative cooling, stellar winds, supernova 
feedback and galaxy merging) and represents the endpoint 
of galaxy evolution along the cosmic history. The next step 
would be the investigation of the evolution of the population 
gradients with redshift, which will be the topic of a forth- 
coming paper. Even at this stage though, we can build a 
clear picture of the relative contributions of each physical 
process in shaping the present day trends. 

The overall trend of metallicity gradients with mass is 
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only partially consistent with the expectation from a mono- 
lithic collapse, where the accretion of SSPs produces a higher 
central metallicity and negative metallicity gradients, that 
are steeper at larger masses where the potential well is 
deeper. 

At high masses, the null age gradients are the conse- 
quence of galaxy merging that produces a mixing of the SSPs 
after the event. The null gradients are in fact more common 
in cluster environment and contrast with the presence of 
positive gradients in (normal) group environment (see e.g. 
Fig. [1} , where there are still a few massive star forming mas- 
sive systems (Fig. [2} . The "merging remnant" systems are 
also the ones having the shallower (negative) metallicity gra- 
dients with respect to the steeper ones of the star-forming 
systems. This shows that at the massive side merging is still 
a major player of the galaxy evolution -what makes our re- 
sults consistent with former simulations of galaxy mergers 
(see Fig. 1 again). 

In fact, typically major dry mergers are known to 
level out pre-exist i ng metallicity and colo ur gradients (see 
iPipino et all l20ld |Pi Matteo et all 120091 ') . hence this fea- 
ture imposes a constraint upon the relevance of merging over 
passive evolution at the high mass end. In particular, at the 
highest masses the occurrence of strong mass accretion due 
to minor and, mainly, major mergings prevents the gradients 
to further steepen. 

At these mass scales, the amount of energy ejected by 
SNe contributes to yield steep Z gradients in massive star- 
forming galaxies, which still have to experience major merg- 
ing events. On the other mass side (i.e. logM, < IOMq) 
instead, SNe produce, on average, null gradients due to the 
large meshing action required by winds easily propagating 
within their smaller volume and less deep potential wells. In 
this case we see that star formation has been shut off in all 
classes of galaxies (Fig. [2]), and the age gradients are simi- 
larly close to zero. Environment acts as to make galaxies in 
clusters have negative gradients with respect to the groups. 
This suggests a different response of the enriched medium 
against SNe explosions. In lower density environments, the 
metals are kept in the small systems and recycled until the 
overall metallicity gradients are swept out, while in the clus- 
ter core environment the higher density and the more nu- 
merous encounters imply that part of the metals are lost in 
the inter-galactic medium and the net metallicity gradients 
stay negative as for massive systems. 

However, there are a number of galaxies (both in clus- 
ter outskirts and groups) which seem to diverge from the 
overall age trend at the low mass end, since they show pos- 
itive age gradients but metallicity gradients consistent with 
the averages of the sample. These systems do not corre- 
late with other galaxy properties, like the star formation, 
color or velocity dispersion (see e.g. Fig. [2]). Instead they 
are the systems with lower ages (see Fig. [3}, well matching 
the suggestion that age is responsible of the scatter of the 
gra dients (especially th e age ones) in the observed systems 
fsee lTortora et aUl2010l and Fig. [3j. These are systems that 
have born in the late epochs and that have evolved quite pas- 
sively so far without experiencing any other stellar processes 
but stellar ageing. We have checked that these galaxies are 
gently migrating toward the zero gradients which they will 
reach in the near future. Most of these dwarf-like systems 
resulted to be star forming at z > (R+08) and are now 



dead mainly as a consequence of their interaction with close 
companions (as demonstrated by the fact that the few star 
forming galaxies are in normal groups and cluster outskirts), 
while SNe seem not to play a major role. At the same time 
these systems could not be residing in FG and inner cluster 
regions because they had merged at earlier epochs onto the 
brightest central galaxies. 



5 CONCLUSIONS 

In this paper, we have used the N-body+hydrodynamical 
simulations described in R+08 to study the age and metal- 
licity gradients in galaxies as functions of mass in the range 
log M* G (9.5 — 11.5) and environment (from groups to clus- 
ter centers). Our results show different trends of age and 
metallicity gradients with mass and galactic environment 
at z—0. With respect to the mass: 1) dwarf galaxies have 
generally null age gradients and shallower or null metal- 
licity gradients (with few cases of positive age gradients); 
2) massive early-types (logM* ~ 10. 5Mq) have moderate 
age gradients that slightly increase with mas s, and nega- 
tive m etallicity gradients (~ —0.2, —0.3, see e.g. lRawle et al.l 
2010, T+10), which reach the steepest values ~ —0.4 in few 
of the most massive galaxies (logM* ;>10.5). At fixed stel- 
lar mass there is a clear dependence on the central galaxy 
age, since younger galaxies have positive age gradients and 
slightly steeper Z gradients, and central metallicity. Simi- 
larly, some marginal trends are found as a function of SSFR 
and colour: in particular the very massive systems with bluer 
colours and active SF are found to have steeper gradients. 
This mass separatio n in the profile slope is in agreement with 



1 nis mass separation in ttic promc slope is m agrccmcn 
other models (e.g. Kawatal l200ll . IKawata fc Gibson! 



Hopkins et al.ll2009a), and observations |Rawle et al.l 



2003, 



2010, 



Ispolaor et"aiTl2009l . iTortora et afll2010l ). even though the 
statistic for our very massive side is quite poor. 

Such behaviour can be explained in terms of the differ- 
ent role played by merging. At log M* < 10.5 mergers are 
rare and thus less effe ctive in mixing stellar population (e.g. 
Ide Lucia et al. Il2006j ): then the only physical mechanisms 
in action are the SN feedback and the passive ageing (simi- 
larly to a monolithic collapse scenario), t hat result into the 
decreasing trend of Z gr adients (see also IKawata fc Gibson! 
120031 . IPipino et al.ll201oh . 

On the other side, the flattening of metallicity gradi- 
ents trend at larger masses can be due to the higher fre- 
quency of major dry mergers at low redshift, which have an 
increased efficiency in producing flatter metallicity profiles. 
IPipino et all |2010) have demonstrated that equal mass dry 
mergers between ellipticals systematically halve the slope of 
any pre-existing metallicity gradient. 

With respect to the environment, systems with posi- 
tive age gradients tend to be found in the external parts of 
clusters and in groups, independently of galaxy 's mass and 
with a quite large spread; while galaxies in the cluster cores 
and fossil groups have, on average, quasi-null age gradients, 
again at all masses. However, when analyzing the average 
values, cluster dwarf galaxies have, steeper metallicity gra- 
dients with respect to the ones in groups, which also present 
shallower values more peaked around zero (with a larger 
fraction of positive gradients); slight steeper age gradients 
are observed in cluster galaxies too. At very high mass, al- 



Stellar population gradients 7 



beit the lower statistics, this trend gets inverted and galax- 
ies in groups have steep er (negative) metallicity gradients 
l|La Barbera et aD 12005^ 1 . As to the scatter of the relations, 
cluster cores and FGs present the tightest distribution of age 
gradients, while the metallicity gradients show no strong dif- 
ferences in the scatter among the different environments. 

Therefore, besides the main role of SN feedback, we have 
found that environment shapes differently the gradients at 
low and high masses. At low masses, tidal interactions are 
acting to give steeper metallicity gradients in cluster galax- 
ies, while merging produces shallower gradients in massive 
systems in the same environment. 

In forthcoming works, we plan to extend this analysis to 
higher z and as well to the BCGs, to put on a firmer ground 
our physical interpretations. In particular, we aim at further 
investigating the physical processes that rule the very mas- 
sive g alaxies such as galaxy major merging and AGN feed - 
back |Antonuccio-Delogu fe slik1l200St iTortora et al.ll2009l ). 
which would be important to shape both the global proper- 
ties and the gradients in stellar populations. 
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APPENDIX A: SYSTEMATICS 

In this section we will analyze various systematics induced 
by selection criteria and the fit procedure adopted to derive 
the age and metallicity gradients. 

Al Fitting procedure 

As described in the text, we have performed a linear fit 
log X — log R in order to recover the age and metallicity pro- 
files of the simulated galaxies. The first basic step of such 
a procedure is to derive the galaxy centre from the distri- 
bution of the star particles. The particle distributions for 
two template galaxies are shown in Fig. IA1I The choice of 
the median as the galaxy centre is dictated by the scarce 
sensitivity of such estimator to the outliers: in many cases 
indeed, the particle distribution is not symmetric (as shown 
in the bottom panel of Fig. IA1|) . therefore the median is 
more suitable to determine the minimum of the potential 
well. 

The profiles for galaxy age and metallicity are shown in 
Fig. IA1I too. The data points are compared with our best 
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Figure Al. Particle distributions, age and metallicity profiles from the left to right, for a massive (top panels) and dwarf (bottom 
panels) group galaxy. Left panels. Distribution of particles in the plane x— y. The axes coordinates are re-normalized to the the median 
of the particles (red point), while the blue point is the mean. The red contours show the density of data-points. Middle panels. Galaxy 
age profile shown as log age vs log R, the contours show the density of data-points, the blue bars are the medians with the 25- 75th 
percentiles, while the red line is the linear fit we have performed. Right panels. The same as the middle panels but for the metallicity 
profile. 



fitted profile, iso-density contours and medians in different 
radial bins. The fitting procedure is directly performed on 
the collection of points, cutting the tails of age and metallic- 
ity distributions out of 3u, to avoid systematics from outliers 
affecting the slope's estimate. 

Finally, we have also analyzed the change in the derived 
gradients when the age and metallicity profiles are fitted in 
a limited radial range, as made with observations. We fit the 
profiles using only the particles in the range (_R o ff/10, R c s), 
where R e g is the radius which encloses half of the stellar 
mass, and the final results are qualitatively unchanged. 

A2 Systematics induced by the small number of 
particles 

We have performed a series of tests on our galaxies to un- 
derstand what is the spurious scatter introduced in fitting 
objects with a very low number of particles. We have taken 
some galaxies with a high number of particles (in the range 
1000 — 5000) and fixed the derived slopes as with a null un- 
certainty. Then we have extracted from them 1000 synthetic 
galaxies having a number of particles N par = 70, 225, 710 
and when possible = 2250, which correspond, for the high- 
est resolution in our simulations, to log M* ~ 9.5, 10, 10.5 
and 11.5, respectively. The fitting procedure is performed 
on each of these 1000 synthetic systems and for each N par ; 
from the slopes hence obtained we have derived a best me- 
dian value and a la uncertainty. While the best fitted slope 



of the original galaxies are perfectly recovered, with an al- 
most null scatter, the uncertainty is not negligible and is, 
obviously, larger for very low mass systems, while is very 
small for the most massive galaxies with logiW, ~ 11. We 
have presented these results in Fig. [T] 

As already discussed, this uncertainty can produce part 
of the larger scatter observed in the gradients at very low 
masses and these systematics could be stronger when the 
number of particles is lower than the minimum value of 80 
we have imposed. To understand the impact on our results, 
in Fig. IA2I we show the age and metallicity gradients for 
both group and cluster members as a function of mass, clas- 
sifying them on the basis of the number of particles. The 
cluster members are the systems strongly affected by our 
selection criterion on N par , that cuts out many galaxies 
with logiW, < 10.5. Such a cut reduces the spread in the 
metallicity gradients, while for the age gradients no relevant 
variations are reported. 
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Figure A2. Age (left panel) and metallicity (right panel) gradients as a function of stellar mass. Cirles and boxes are for cluster and 
group members, respectively. Filled and open symbols are for galaxies with N par > 80 and N par ^ 80, respectively. 



